ETH-TH/94-20, |hep-th/(J411I^ 



> 



Oh: 



s 



Normalization of scattering states and Levinson's 
theorem: reply to a comment by 
[Helv. Phys. Acta 67, 20 (1994)] 



o\'. theorem: reply to a comment by R. G. Newton 



> 

o 

t^ ; By Nathan Poliatzky 



Institut fiir Theoretische Physik, ETH-Honggerberg, CH-8093 Zurich, Switzerland^ 



Abstract. Normalization of scattering states, Levinson's theorem and the high energy behavior of 
scattering phase shifts are discussed in the light of R. G. Newton's recent criticism. 



j>^ '. 1 Introduction 



Recently R. G. Newton published a comment 0] criticizing the methods and the results of a 
paper published in this journal 0. His criticism touches on a few key points of the subject 
and hence deserves a detailed reply. Here is the reply, point by point, to his criticism. 



2 Normalization of scattering states 

Let us first consider the case without a threshold resonance (i. e. either a finite solution 
with k = does not exist or it is a bound state). We will return below to the (rather 
exceptional) case with a threshold resonance (sometimes called a half-bound state or a zero- 
energy solution). In this section we consider the case of a Schrodinger equation; the case of 
a Dirac equation is completely analogous and the relevant formulae can be found in 0. 
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The usual statement on the normahzation of scattering states for a central potential is 

dr Uki {r) Uk'i (r) = 27t6 {k - k') (2.1) 



/o 

The question then arises of what happens in the limit k -^ k' . Of course the right-hand 
side becomes infinite, but it is easy to see that the difference, between this equation and its 
free-particle counterpart, is finite. In [^ it was proven that 



oo 

2 r^\ „.2 



dr 





nil (r) - vl (r)J = 2r^[ (k) (2.2) 

where Vki (r) are the free-particle wave functions and rji {k) is the phase shift. 
Equations ( p.l| ) and ( |2.2| ) can be put together as one compact equation 

("OO 

/ dr Uki (r) Uk'i (r) = 2n6 {k - k') + 2r]l {k) 6kk' (2.3) 

Jo 

where 6kk' = 1 ii k = k' and 6kk' = otherwise. 



Equation ( ^.31 ) is what R. G. Newton criticizes. As he put it: "He arrives at a modifi- 
cation of the usual quasi-orthogonality of the scattering wave functions in which the Dirac 
distribution is supplemented by a bounded function that differs from zero at only one point. 
In any context in which the Dirac distribution is meaningful, such a function is of course 
equal to zero, and its addition is of no consequence." 



Admittedly, if the delta function on the right-hand side of ( p.3|) is interpreted strictly 
as a distribution, the addition of the second term is of no consequence and equations { \2.'S\) 
and (|2.1| ) are identical. However, the above criticism misses several points. One is that ( p.3|) 
is merely a short-hand description of the two equations (|2.1D and (|2.2D, the validity of which 



is beyond any doubt. So that, if needed, any reference to ( |2.3|) can be avoided. 



A deeper point is that the Dirac delta function and a distribution are not exactly identical 
concepts, e. g. the delta function on the right-hand side of ( p.l| ) or (|2.3| ) cannot be interpreted 
as a distribution in the usual mathematical sense. The proof of this statement is as follows. 
Assume that this delta function is a distribution. However, for the left-hand side of ( p.l|) 
it perfectly makes sense to subtract its free-particle counterpart and then put k = k' . The 
result is finite and given in ( p.2|) . On the other hand, for the right-hand side of ( p.l|) such 
an operation does not make sense (or its outcome is ambiguous), if the delta function is 
interpreted as a distribution. Hence the left- and the right-hand side of ( |2.1| ) cannot be the 
same object and the conclusion is: the right-hand side of equation ( p.3|) is a modification of 
the usual notion of a distribution. 

The mathematical notion of a distribution is based on a class of all possible limiting 
processes, whereas in quantum theory the Dirac delta function is always either a particular 
limiting process or a particular class of limiting processes with a fixed physical interpretation 
of the parameters involved.^ Having this in mind there is no danger of ambiguities for the 



use of equation (|2.3| ) 



^For instance, in the usual calculation of a differential cross section in quantum field theory one has to 
square a Dirac delta function, which is done unambiguously using a particular physically motivated limiting 
process. From the point of view of a distribution this does not make sense. 



3 Levinson's theorem 

Consider the completeness relation 

POO Aj^ 

V u,^i (r) u,,i (r') + / --Uki (r) Uki (r') = 5 {r - r') (3.1) 

where the sum runs over the bound and the integral over the scattering states. 



If the right-hand side of (|3.1|) is interpreted as a Dirac delta function (i. e. as a class of 
limiting processes which are independent of the potential) and not as a usual distribution, 
we can subtract from ( p.l|) its free particle counterpart and put r = r'. The delta functions 
cancel and we obtain 



E <i (0 + / TT- [nil (r) - vl (r) = (3.2) 

^n JO ZTT L J 



dk 
Ki [r) + / 

E,<0 

Integrating (|3.2| ) over r and substituting (2^), we obtain the Levinson theorem 



VI (0) = ri^TT (3.3) 

for the case of a Schrodinger equation, where ni is the number of bound states. If the same 
is carried out for the case of a Dirac equation (see ref. [@] for the details), the result is 

ri^. (0) + r/_„,. (0) = (iV+ + N-) ix (3.4) 

where N'^ is the number of positive and A^~ is the number of negative energy bound states.^ 



Looking at ( |3.4|) it is natural to ask whether a stronger statement of Levinson's theorem is 
possible, valid for the positive and negative energy phase shifts separately. It turns out to be 
the case Q. The clue to derive such a stronger statement of Levinson's theorem is simple: 
consider the second order (iterated) Dirac equation and take the limit /c — ;> 0. In this 
nonrelativistic limitQ the second order Dirac equation becomes identical to a Schrodinger 
equation, so that at A; = one can relate the phase shifts of a Dirac equation to the 
phase shifts of a corresponding Schrodinger equation. The stronger statement of Levinson's 
theorem obtained in this way is: 

VniK (0) = n+TT (3.5) 

?7_„«, (0) = n^TT (3.6) 

■^More precisely we obtain 

m (0) = -qi (oo) + niTT 

for the case of a Schrodinger equation, and 

r/mK (0) + Tl-,n^^ (0) = 7700,. (00) + 77_oo.«: (oo) + {N+ + A^^") TT 

for the case of a Dirac equation. However, we are free to put 77/ (00) — and tjook (00) + ??-oo,k (oo) = 0. 

*It is an interesting fact that such a non-relativistic hmit which, from the physical point of view, is a most 
natural one, does not coincide with the usual expansion based on the Foldy-Wouthuysen scheme (see B). 



where ri^ is the number of positive and n~ the number of negative energy nodes of the A; = 
solution of the Dirac equation. Unhke the case of a Schrodinger equation, there is a subtlety 
here: for a given energy sign the number of bound states may not equal the number of nodes 
of a /c = solution, i. e. in general A^^ 7^ ri^ and A^~ 7^ n~. The equality holds only for the 
total numbers: 

Nt + N-=nt + nl (3.7) 

The generality and simplicity of the above proof of Levinson's theorem speaks for itself, 
and can be viewed as a confirmation of the validity and usefulness of the above discussion 
on the normalization of scattering states. Notice that we did not have to assume anything 
explicit about the potential. The only implicit assumption is that the potential decays faster 
than 1/r^ at large distances, which is sufficient for the validity of (|2.2|) . 



However, when in ref. [@] this proof was extended to the case with a threshold resonance, 
a surprising result was obtained: 

m (0) = {ni + g) TT (3.8) 

in the case of a Schrodinger equation, and 

V±mn (0) = {n^ + g) vr (3.9) 

in the case of a Dirac equation, where either g = 0, org = i, org = |. For the case with a 
threshold resonance g = can be ruled out and g = | is what one would have expected. The 
surprise is in the possibility of allowing for q = \, which is the subject of R. G. Newton's 
criticism. From the point of view of a conventional proof of Levinson's theorem, where 
there is always a certain restriction on the potential, the case q = \ cannot be ruled out 
completely. However while the paper was in press it was found that the derivation of the 
normalization integral actually implies a certain constraint which was overlooked in and 
which rules out the case q = \- This is the reason why in a subsequent paper the case 
q = \ was not mentioned. Meanwhile M. Sassoli de Bianchi investigated the problem 
within the framework of one spatial dimension and came to the same conclusion. In the 
following section these matters are considered in more detail. 



4 The case with a threshold resonance 

Assume that the coupling constant is tuned such that a finite, fc = solution (threshold 
resonance) exists. By examining its asymptotic behavior at large distances, it is easily seen 
that such a solution is normalizable (i.e. a bound state), unless / = 0, in the case of a 
Schrodinger equation, and k = —1, e = m or k = 1, e = —m, in the case of a Dirac 
equation.0 Hence in this section / = 0, in the case of a Schrodinger equation, and k = ±1, 
in the case of a Dirac equation. 

^Notice that, for a given angular momentum, at most one threshold resonance is possible, either 
with e = rn or with e == —m. 



According to [^ in the presence of a threshold resonance ( |2.2| ) becomes modified according 
to 



dr 





«L (r) - vl, (r) = 2i, (k) + 2n6 (k) sin^ r/o (A;) (4.1) 



As in section ^, integrating ( |3.2| ) over r and substituting ( [4.1| ), we obtain Levinson's 



theorem for the case with a threshold resonance 

770 (0) = noTT + I sin^ 7^0 (0) (4.2) 

Equations ([4.1| ) and ( |4.2| ) are valid for a Schrodinger equation. In the case of a Dirac 
equation, as was mentioned in the preceding section, one can rewrite the Dirac equation at 
/c = as an ordinary Schrodinger equation and use ( |4.2| ) to obtain the stronger statement of 
Levinson's theorem for a Dirac equation: 

V±mK (0) = n^rc + - sin^ ?7±„« (0) (4.3) 

where ra^ is the number of positive and n~ the number of negative energy nodes of the k = 
solution. 



Equation ( [4.2|) was first derived in ref. 0. The authors did not check for all solutions 



of this equation but instead were content by verifying that (|3.8| ) for q = and g = ^ is a 
solution of ( ^.2| ). Equations ( |4.2| ) and ( [4.3| ) were derived independently in ref. and it was 



pointed out that there is a third solution with q = \- However it was not noticed that the 
derivation of ( |4.1| ) also implies 

sin [2r]i (0)] = (4.4) 

This constraint rules out the case g = | in (^]8|).0 



To understand how this constraint comes about, we now consider the derivation of ( [4.1D 
in more detail (section 2.1 of ref. 0). We start with the exact identity 

l"" dr uli (r) = ^ [uii (R) duUu (R) - uui (R) duu',, {R)] (4.5) 



and assume that the potential V (r) vanishes for r > R. Then the right-hand side of ( [4.5|) 
depends only on the asymptotic wave functions and thus can be calculated. The result is 

f dr uli (r) = 2R + 2r][ {k) - (-1)' - sin [2kR + 2rii (k)] (4.6) 

Subtracting from (|4.6|) its free-particle counterpart and expanding the sinus, we obtain 



^The situation for a Dirac equation is completely analogous and the constraint is 

sin [2r]±^^ (0)] = 



This constraint rules out the case q — j in (3.9) 



R 



dr 



uli {r) - vli [r] 



H ik) + i-iy 1 2 sin^ VI {k) ^^^^ - "" ^?' ^^^^ cos {2kR) ' 



k 



k 



(4.7) 
The second term in curly parenthesis possesses a k ^ singularity unless sin [2rii (0)] = 0. In 
the case with a threshold resonance, where the k = state is part of scattering states, such a 
singularity is intolerable (for instance, the integral over k in ( p.2| ) includes the point k = in 
this case). Hence the constraint ( [4.4|) must be satisfied. The derivation of ( ^4.1| ) is complete 
now, if one sends i? — »• cxo. 



5 High energy limit 



In it was shown that the high energy behavior of the scattering phase shifts is 



Vi{k) 



fc— »oo 



oo 

* k I dr 






(5.1) 



in the case of a Schrodinger equation, and 



VeK {k) 



fc^oo 



k I dr 



'i-|vw-i 



(5.2) 



in the case of a Dirac equation. Obviously, if the potential V is less singular than 1/r at the 
origin, then (|5.1|) and ( ^.2|) lead to the well known results 



7]i (oo) = 
in the case of a Schrodinger equation, and 






(5.3) 



(5.4) 



in the case of a Dirac equation. 

R. G. Newton doubts the validity of (|5.1|) and ( ^.21) : "His change of variables leads, 
without comment, to a highly singular equation on a complex contour, and it is not clear 
whether his manipulations are valid." 



This criticism is without foundation, since the equations, which were used to derive ( p.l|) 
and ( |5.2| ), become singular only at certain finite values of k. Taking into account the fact 
that these equations are local in fc, such singular behavior is irrelevant to the high energy 
limit k ^ oo. 
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